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Abstract — Pulsatile flow of a Jeffrey fluid In a channel 
bounded by porous lined plates with suction and Injection 
is studied in this paper. The steady and unsteady 
velocities are obtained. The effect of various parameters 
on the flow phenomenon is discussed through graphs.. 
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I. INTRODUCTION 

Viscous fluid flow through and past porous media is 
attracting the attention of scientists and engineers 
because of its wider applications in various branches of 
science and technology. The movement of ground water 
in soil, the seepage of water through earth fills and 
concrete dams, the movement of oil fields can described 
using the knowledge of flow through porous media. The 
petroleum industry has been showing a lot of interest in 
these problems in connection with the crude oil 
production from the under ground reservoirs. These 
reservoirs consists of porous materials like lime stone 
and dolomite where oil is preserved. Oil can be obtained 
by drilling wells down into the reservoir. In order to have 
a better oil production, it is necessary to use the 
knowledge of flow through porous media. 



flows through and past permeable beds. 
Radhakrishnamacharya [6] investigated the pulsatile 
flow of a dusty fluid containing small solid particles 
uniformly, through a two-dimensional constricted 
channel. The effect of non-Newtonion nature of blood 
and pulsatality on flow through a stenosed tube is 
analyzed by Chaturani and Samy [7]. 

The pulsatile flow in a porous channel is important in 
understanding the process of dialysis of blood in an 
artificial kidney. Recently, Chandra and Prasad [8] 
discussed the pulsatile flow problems with periodic 
acceleration and varying cross section of tubes. 

Wang [9] studied the interesting problem of pulsatile flow 
in a porous channel bounded by rigid walls. The 
pulsatile flow between permeable walls is important in 
understanding the blood flow in the circulatory system 
where the nutrients are supplied to tissues of various 
organs and waste products are removed. Vajravelu et 
al.,[10]made a detailed study on pulsatile flow between 
permeable beds. Avinash et.al [11] studied the pulsatile 
flow of a viscous stratified fluid of variable viscosity 
between permeable beds is studied. The interaction of 
peristaltic flow with pulsatile flow through a porous 
medium is discussed by Afifi and Gad [12]. 



The oil available in the porous reservoir is a complex 
fluid. The properties of such fluid have impact on the oil 
production. The behavior of the oil may be Newtonian or 
non-Newtonian. In view of this, it is interesting to study 
non- Newtonian fluid flow through and past porous 
media. Further there are many important applications in 
biomechanics also (vide Fung and Tang, [1,2]. 

Muakat [3] made theoretical and experimental studies 
on porous flow using Darcy law. Darcy law is observed 
to be valid for low speed flows and agrees with several 
experiments modeled in one dimensional motion. Yih [4] 
suggested the modified Darcy law for describing 
unsteady flow through porous media. Following this law, 
Rudraiah et al. [5], analyzed several time dependent 



n this paper an exact solution for the pulsatile flow of 
Jeffrey fluid in porous lined plates is obtained. The flow 
between the permeable layers is governed by Jeffrey 
model where as the flow in the lower and upper beds 
are governed by Darcy law. The velocity distributions in 
the porous and non-porous regions are determined. 

II. MATHEMATICAL FORMULATION OF THE 

PROBLEM 

Consider the pulsatile flow of a Jeffrey fluid in a 
channel bounded by porous lined plates (Fig. 1 ). 
The thickness of the porous lining on each of the 
plates is e. The fluid is injected into the channel 
from the lower porous layer with a velocity V and is 
sucked out into the upper porous layer with the 
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same velocity. The permeabilities of lower and 
upper beds are and k 2 . The flow between the 
permeable beds is governed by the Jeffrey model 
whereas the flow in the porous medium is 
described by modified Darcy's law 



du 



a 



(Mbi ~Qi) 



u — u 



B2 



du -a 



dy Jk~ 2 



(u B2 -Q 2 ) 



at y = E 

at y = h - E 
at y = h - E 



(5) 

(6) 
(7) 



h-2 E I 



Centre Line 



u, v are velocities of the fluid, v is the coefficient of 
viscosity 

p is the pressure, a is the slip parameter, k 1 ,k 2a are 
permeabilities of the lower and upper permeable 

beds, u Bl9 u 52 are slip velocities in lower and upper 
permeable beds and \ is the Jeffrey parameter. 



Fig. 1 . Physical Model 



The following assumptions are made in the 
analysis of the problem 

The flow is laminar and fully developed 

The permeable beds are homogeneous 

The flow is driven by unsteady pressure gradient. 
We assume thatj_3p =>4+& to 

p dx 

where A and B are constants and 1 a)' is the 
frequency. In view of the above assumptions, the 
basic equations and boundary conditions of the 
flow take the following form. 



Basic equations 



du dv 
+ — = 0 



dx dy 



du Tr du -I dp v d 2 u 
— + V — = — + - 

dt 



(1) 



dy p dx \ + \ dy' 



^ = 0 

dy 



(2) 



(3) 



Boundary conditions 



u — u 



Bl 



at y = E 



(4) 



Separating equations (1) - (7) into steady 
part denoted by a bar (-) and unsteady part 
denoted by a tilde (~) we get 



Steady part: 

du =Q 



dx 



u d u 



V — = -A + 
dy \ + \ dy 2 

u — Ub\ 



du a - — 
(ubi -&) 



dy 



u = u m 



aty = E 



aty = E 



at y = h - E 

^=" (UB2 -Q 2) aty = h-E 
dy ^jk 2 

where q q h^L^ K ) 



Unsteady part: 

^ = 0 

dx 

du | du _ Beiat | v d 2 u 
dt 
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du -a 



dy Jk~ 2 



(u B2 Q 2 ) at y = h - e 



(19) 



where g = -k lP B{\ + \) ^ > ^ = -k 2 pB(l + ^) ^ 

Non - Dimensionalization of flow quantities 

The following non - dimensional quantities 
are introduced to make the basic equations and the 
boundary conditions dimensionless 

Steady part: 



u 



u = 



(4*) 

V 



— * Ubi 

Ub\ = 



4* ■ 



— * UB2 
UB2 = 



4* 



X — — j 

h 



V 



V 



V 

s = — where X = . 
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Unsteady part: 



Q 



2 » 



V 



U 
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- 



51 



1/ 
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52 
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r~ — ' 

e, 
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■ — ^ 
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, y = — , £ = — where 5, = -5 . 

A A 1 



1/ = 1/52 

du 



at y = 1 - 8 



dy 



1 _ i - ,/L 

= -«o- 2 (wij2 — — ^-i?) aty = 1-£ 



(24) 
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where = — is the Reynolds number, <j 1 = 



v 



(To = 



(o" 19 o" 2 are dimensionless parameters). 



Unsteady Part: 

du 



dx 
d 2 u 



= 0 



-£(1 + 4) — -(1 + 4) — = -(l + 4> 



Letting w = f(y)e ia * 

f\y) - R{\ + 4 )/' (y) - / 0(1 + A )7(y) = -(1 + A ) 

/— — ' 

Using w sl = j x e , u B2 = J 2 e 
The boundary conditions become 

f = fi at y = £ 

J/ - ~ 1 + A, . 
-j- = aa l (f l ^-) at y = £ 

dy 

f = f 2 

df 



(26) 



(27) 



(28) 



at y = 1 - £ 



= -aa 2 (7 2 -^ 1 ) aty = 1 - 
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I. SOLUTION OF THE PROBLEM 



In view of the above dimensionless quantities 
Eqs.(8) to (19) take the following form, neglecting 
the asterisks (*), we get, 



Steady part 



Steady part: 



Solving Eq.(21) subject to boundary 
conditions (22)-(25) we get the velocity field as 



We get the velocity field as 

u = C l +C 2 e (l+Ai)y + y 



(33) 



du 



dx 
d 2 u 
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= 0 



du 
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dy 

U = Ubi at y = e 

du - 1 + 4 
— = ao x (u b\ r—7?) 
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at y = £ 



(20) 
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(22) 
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The slip velocities are Ubi and u B i are given by 

- {D 2 E 2 -D,E X ) - _ {D,E X -D X E 2 ) 

Ub2 — ~ (34) 



Ubi = 



(AA-AA) 



(AA-AA) 



Unsteady part: 



Solving equation (28) subject to boundary 
conditions (10) -(13), we get 
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f(y) = Ae m ^ y +Be m2y - 



CO 



(35) 



The unsteady part of the velocity is given 

u = f(y)e iwt 

Separating real and imaginary parts, we get 

u — (A l 2 Cos cot — B x 2 Sin cot) + i(B l 2 Cos cot — A 12 Sin cot) 
The slip velocities are given by 

^ 1 + 4 X i 

-+ — T-y 



by 



u 
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5. Deductions 

(i) Taking k 1 = k 2 = k (i.e. a 1 = o 2 = o) in equation (33 ) 
and (35) we obtain the velocity field for pulsatile flow of 
Jeffrey fluid between porous lined plates with equal 
permeability as follows: 



Steady state velocity 



u- + c 2 e 



(l+X)Ry 



(36) 

Unsteady part: 

u = (A 1 e m ^ + A 2 e m ^ 

(37) 



M 2 J 



icot 



(ii) When the permeabilities ki and k 2 tend to zero 
in equations (36) and (37) we obtain the velocity 
field for the pulsatile flow of the Jeffrey fluid in a 
channel bounded by rigid walls as follows 



Steady part 

u = C ± + C 2 e Ry + y 

The slip velocities u B1 and u B2 are zero. 
Unsteady part 

u = (A^y + A 2 e m *y - j^j e ia)t 



(38) 



(39) 



The slip velocities u B1 ar\d u B2 are zero. 

Further with A -> 0, the results (38) and (39) reduce the 

corresponding ones of Wang [9]. 



RESULTS AND DISCUSSION 

From equation (33) we have 
calculated the steady part of the velocity as a 
function of y for different values of permeability 
parameter o with fixed a = 0.5, R = 2, A = 0.5, 

£ = 0.02 and is shown in Fig. 2. We observe that the 



velocity decreases with the increase in G = — . This 

may be due to the increase in the permeability 'k' of 
the upper and lower beds. Also, as the a 
increases, the gap between the velocity profiles 
becomes smaller i.e. there is not much change in 
the velocity due to variation for large values of a . 

The variation of steady part of the 
velocity with y is calculated from equation (33) for 
different values of permeability parameter a with 

fixed a = 0. 1, R = 2, X = 0.5 , £- 0.02 and is shown 

in Fig. 3. It is noticed that the velocity decreases 
with the increase in a . From Fig. 2 and Fig. 3 it is 
noticed that the magnitude of the velocity increases 
with decrease in the slip parameters . Also, as the 
a increases, the gap between the velocity profiles 
becomes smaller. 

The variation of u with y is 
calculated from equation (33) for different values of 
the Reynolds number R and for fixed 
a = 2, cj = 25, /L = 0.5, s = 0.02 and is shown in 

Fig. 4. We observe that the velocity increases with 
the increase in R. 

From equation (33) we have calculated the 
steady part of the velocity as a function of y for 
different values of Jeffrey parameter a with fixed 

a = 2, cr = 25, R = 2, £ = 0.02 and is shown in 

Fig. 5. It is noticed that the velocity increases with 
the increase in the Jeffrey parameters. Similar 
behavior due to a is noticed by Srinivas et al. 
(2008) and Hayat et al. (2008) for the peristaltic 
transport of viscous fluid in a flexible channel. 

From equation (33) we have calculated the 
steady part of the velocity as a function of y for 
different values of the thickness of porous lined 
plate s with fixed a = 2, a = 25, r = 2 ,1 = 0.5 and is 

shown in Fig. 6. We observe that the velocity 
decreases with the increase in s . 

We have calculated the unsteady part of the 
velocity as the function of y from equation (35) for 
different values of wt with fixed a = 5, cr = 5, R = 0 

,^ = 0.02, M=1, yl = l and is shown in Fig.7.We 
observe the velocity decreases with the increase in 
wt . 

We have calculated the unsteady part of the 
velocity (u ) as the function of y from equation (35) 
for different values of wt with fixed 
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a = 0.5, cr = 25, /? = 10,<? = 0.02, M = 1, A = \ and is 

shown in Fig.8.We observe the velocity decreases 
with the increase in wt . From Fig. 7. and Fig. 8. for 
fixed wt and a, as R increases the maximum 
velocity moves closure to the upper permeable 
bed. 



From equation (35) we have calculated the 
unsteady part of the velocity as a function of y for 
different values of Reynolds number R and Jeffrey 

parameter a with fixed a = 5, M=1, s = 0 02 wt = — 

4 

, cr = 10and are shown in Fig. 1 1 . and Fig. 12. It is 
noticed that the velocity increases with the increase 
in a . As R increases the maximum velocity moves 
closure to the upper permeable bed. 

The variation of a with y is calculated from 
equation (35) for different values of the Reynolds 
number R and for fixed a = 2, <y = 25, M=1 , 

z = 0.02, = — , k = 0.5 and is shown in Fig. 13. We 

4 

observe that the velocity decreases with the 
increase in R. 

From equation (35) we have calculated the 
unsteady part of the velocity as a function of y for 
different values of the thickness of porous lined 

plate 8 with fixed <^ = 2, cr = 10 , M = 1 , R = 2wt = E 

4 

, A = 0.5 and is shown in Fig. 14. We 

observe that the velocity decreases with the 
increase in s . 
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Fig. 2. Steady state velocity profiles for 

a = 0.5, R = 2, A = 0.5 and £ = 0.02 
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Fig. 3. Steady state velocity profiles for a with 
fixed a = 0. 1, R = 2, X = 0.5 and £ = 0.02 




y 

Fig.. 4. Steady state velocity profiles for R with fixed 

a = 2, cr = 25, A = 0.5 and £ = 0.02 



We have calculated the unsteady part of the u 
velocity ( a ) as the function of y from equation (35) 
and are shown in Figs. 9 and 10. For fixed wt , the 
velocity decreases with the increment in a. As R 
increases the maximum velocity moves closure to 
the upper permeable bed. 
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Fig.. 5. Steady state velocity profiles for X with fixed 

a = 2, cr = 25, R = 2 and s = 0.02 
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Fig. 6. Steady state velocity profiles for s with 
fixed a = 2, a = 25, R = 2 and X = 0.5 



0 5 
01 


! ^^-"^^ 

■ — ~ ** ~ 


0 


: 


0.4 fit o> 


10 


i 








-01 






✓ 










-0 
















-0 5 

















Fig. 7. Unsteady state velocity profiles for cot with 
fixed a = 5, c>=5, R = 0, £ = 0.02 , M=1 and X = 1 
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Fig. 8. Unsteady state velocity profiles for cot with 
fixed = 0.5, cr = 25, 7? = 10 , £ = 0.02 , M=1 and 
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Fig. 9. Unsteady state velocity profiles for a with 

2 = 0.5 
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Fig. 10. Unsteady state velocity profiles for a with 

2 = 1 
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Fig. 1 1 . Unsteady state velocity profiles for /I with Fig. 13. Unsteady state velocity profiles for R with 



71 



fixed cr = 5, R = 2 ,M=1 , ^ = 0.02, wt = — 

4 

<7 = 10 



and 



7Z" 



fixed a = 2, cr = 25 ,M=1, ^ = 0.02, = — and 

4 

2 = 0.5. 



> 06 



u >04 



I = l 5 



0.4 



'.'.0 



Vj 




1.0 



u 






£ =0 02 


e = 


0 03 X 




J* 


e=0.04 


J 





Fig. 12. Unsteady state velocity profiles for X with Rg 14 Unste ady state velocity profiles for s with 
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fixed ^ = 5, R = 10, M=1 , ^ = 0.02, w/ = — and 
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fixed = 2, cr = 10, M=1 , R = 2,wt = — and 

4 

2 = 0.5. 
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